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In the context of Modified Newtonian Dynamics (MOND) , we study how perturbation of a spher- 
ically symmetric system would affect the gravitational field. In particular, we study systems of 
perturbed and unperturbed spherical shells. For a single perturbed shell, we show that the field 
inside the shell is much smaller than what would be expected from a naive scaling formula. The 
strength of the perturbation field within the shell is screened by the spherically symmetric com- 
ponent of the mass, and is reduced as the spherically symmetric component is increased. For a 
two-shell system, we again show that the perturbed field is screened by the shells, no matter which 
shell's mass distribution is perturbed. The field within the inner shell is most suppressed when the 
inner and outer shells coincide. However, for a very light inner shell, the perturbation to the field 
can be enhanced. The enhancement is typically larger for smaller inner shells, and the perturbed 
field can be amplified by almost a factor of 2. The relevance to the effect of external fields on galaxy 
dynamics is discussed. 
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PACS numbers: 98.80.-k 04.50.Kd 98.65.Cw 95.35.+d 



INTRODUCTION 

Most of the astrophysical computation concerns a lo- 
cal system embedded in much larger background system. 
For example, Solar system is embedded in Milky Way a 
galaxy in a cluster, a cluster to the supercluster etc. In 
GR or Newtonian gravity, the external field effect on the 
local system due to much larger background mass den- 
sity are often neglected based on two simplifying assump- 
tions: that the background mass distribution is spheri- 
cally symmetric and/or the background distribution is 
located very far from the local system. Under the first 
assumption, one would invoke Birkhoff's theorem that 
the field inside a spherical cavity should vanish. (Inclu- 
sion of the local mass distribution may spoil the spherical 
symmetry of the whole system, and in GR, the back- 
ground mass distribution may exert force on the local 
system. However the external field induced in this way 
is negligibly small. See Under the second assump- 
tion, the inverse-square law will suppress the contribu- 
tion of the portion of the background mass distribution 
that is far from the local system. Therefore, even when 
the background mass distribution is only approximately 
spherically symmetric, the external field effect in GR or 
Newtonian gravity would not significantly alter the local 
dynamics. 

It might be naively expected that the external field ef- 
fect is more significant in MOND theory. A small pertur- 
bation of the background mass distribution from spher- 
ical symmetry might be expected to alter the internal 
field configuration significantly. This is because the field 
of a bounded mass distribution scales as ~ - in MOND, 

r ' 

and the field due to the aspherical part may therefore 
persist to greater distances. Of course, the modified Pos- 



sion equation in MOND is non-linear, and it does not 
embody the superposition principle. One cannot treat 
the contribution from the spherically symmetric part of 
the mass distribution and the aspherical part separately. 
The total mass distribution has to be considered in order 
to infer the size of any external field effect. The external 
field effect in MOND has been receiving growing interest 
[T5L Il6j |. It is argued that MOND theory is consistent 
with the observed escape velocity in the solar neighbor- 
hood if the Milky Way is embedded in a constant external 
field of ~ O.Olao- (Here ao is the universal critical accel- 
eration that characterizes MOND, ao — 1.2 x 10 _8 cm/s 2 . 
On the other hand, Wu et al [17( has showed that an ex- 
ternal field of approximately 0.03ao is needed to bind the 
Large Magellanic Cloud to the Milky Way. Since the ex- 
ternal field effect on the local escape speed presents a 
strict test of MOND theory it is of interest to under- 
stand when, if at all, the aspherical perturbation of very 
distant backgrounds could be neglected. 

In an earlier work [l| , we studied the fields inside spher- 
ical shells as probed by non-negligible test masses. We 
showed that when those test masses are placed off-center, 
thereby breaking the spherical symmetry, the accelera- 
tion of the test mass can become comparable to the (ex- 
ternal) surface gravity of the shell. In what follows, we 
consider simplified systems of spherical shells with per- 
turbations in order to better understand how the per- 
turbations of spherically symmetric mass distributions 
translate into perturbed gravitational fields. 



MOND 

MOND 0, 01, 0] is an alternative to the canonical dark- 
matter scenario in which the gravitational force law (or 
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Newton's second law) is modified to account for the miss- 
ing gravitational field. In 1984, Bekenstein and Milgrom 
introduced the Lagrangian formulation of Modified New- 
tonian dynamics (MOND) 10]. The field equation of 
MOND is derived from the Lagrangian 



L = — d 



pip 



(1) 



where ip is the gravitational potential. J-{y 2 ), with y = 
| V^|/ao, is an arbitrary universal function, that together 
with ao, the characteristic scale of MOND, specifies the 
theory. Varing L with respect to ip yields a modified 
Possion equation: 



V • [,i(|W>|/a )VV] = 47rGp(|r|) 



(2) 



where fx(y) = J-'(y 2 ). fx(y) must approach 1 as \y\ ^> 1 
and |y| as \y\ <C 1, in order that the field scale as 
near a spherical mass distribution (the usual Newtonian 
result) and as - far from the mass distribution to explain 
flat galaxy rotation curves. 

One commonly used form of \x is 



\y\ 



(3) 



(However see [8| for different form of p function.) The 
value of do is then given by phenomenological fit. We will 
adapt the value derived by Begeman et al Q in the study 
of external galaxies with high quality rotation curves. 



a = 1.2 x l(T lu m/s 



(4) 



For a bounded mass distribution of total mass M, we 
define a transition radius 



R t = y/GM] 



Of) 
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R t indicates a point at which the Newtonian field approx- 
imately equals ao, and this is about the point at which 
the field switches from Newtonian 1/r 2 to MOND's 1/r. 

We now define three quatntities, MOND acceleration 
gM , Newtonian acceleration g^ and 'naive' MOND ac- 
celeration g s . MOND acceleration comes from the full 
solution of the modified Poisson equation [2j and is de- 
fined to be 



Sm 



-Vip 



(6) 



Newtonian potential <f> is the soloution of a ordinary Pois- 
son equation 

V ■ VcpN — 47rGp(|r|), (7) 
and Newtonian acceleration gw is defined to be 

gw = -V<j> N . (8) 



From the form of the modified Poisson equation, 'naive' 
MOND acceleration g s and Newtonian acceleration gyv 
for a given density distribution can be related by the 
algebraic relation 



K\gs\/ao)gs = Kn 



(9) 



This is the original formation of MOND and sometimes 
incorporated to calculate the galactic rotation curves @, 
B EH • Inverting the relation we will find 



KlgJvIM)) 



(10) 



where v[x) — I (x)x and I{x) = xfi(x). Usually g s is 
not equal to gM since g s in general is not curl-free. (gM 
by defnition has to be curl-free since it is gradiant of a 
potential). Important exception is if gj\r only depends 
on one parameter, which is the case for spherical, cylin- 
drical or plane symmetry, then gM = g s ||. Then if 
the perturbation is small, we might expect that MOND 
acceleration can well be approximated by g s . 



CODE 

The numerical code used in this paper is developed by 
Milgrom [ll[. It is the same code we used in our earlier 
paper [l[ , and it is the implementation of Milgrom's code 
on a spherical lattice. We set the number of angular glid 
N and that of radial glid L to be 120 and 400 respectively. 
N = 120 tlanslates into the angular resolution of 1.5°. 

PERTURB ATIVE CALCULATION 

We first consider a small aspherical perturbation on 
the shell as in FIG. [TJ The center of the shell is taken 
to be the conter of the spherical coordinate for the sake 
of simplicity. The radius of outer shell and its mass is 
fixed to R = OMMpc and M sh = 2.0 x 1O 13 M which 
would be a typical parameters for a gas- rich cluster 0]. 
Let the density distribution of the perturbation and the 
spherical shell be p\ and p s h, respectively, then we will 
write 



Psh = <JshS(r - R) 
Pl = ea sh Y(9)S(r — R), 



(11) 



where <j s h is surface mass density for the shell, and Y(9) 
is some angular function. The unperturbed field ipo sat- 
isfies the modified Poisson equation: 



V ■ M|VV>o|/ao)ViM = ^G Psh 



(12) 



The exact solution for VV>o can be found by applying 
Gauss's theorem. In terms of the quantity u 



GM S , 



W>0 




u 4 + u z r 



(13) 
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Perturbations 
on the Shell 




FIG. 1: This figure shows the mass distribution of the system. 
A shell with M ah = 2.0 x 10 13 M Q and R = Q.&AMpc is placed 
concentric with the origin of the coordinate . A massless ob- 
server is located at r = 0.05-R with an angular displacement 
of 45° from z-axis. An azimuthally symmetric perturbation 
is placed on the shell. For the multiple shell system, a sec- 
ond shell with mass = M2 with radius — R2 will be placed 
concentric with the first shell. 



1 , 1 
~— a u 2 (1 + -u + ...)r . 

for exterior of the shell, and 

Wo = (14) 
for the interior of the shell. The expansion in small u is 



valid for r S> 



GM k 
an 



The solution ip will satisfies the 



following jump conditions at r=R: 



f\Vip oll 



I Win I 



r\r=R 



by aspherical perturbation, and one may assume that the 
'simple' form of \x function is applicable, (i.e. (j,(x) ~ x 
as is the case for the system in the deep MOND regime). 
Then ipi has to satisfy the following equations interior 
and exterior to the shell. 



\ a o / \ «o 

= 4TrG[a sh + ea sh Y(6)} 

VVW|| - Vlpin\\ = ( 15 ) 

Exterior to the shell, the spherically symmetric solution 
VV'o w ill dominate, and one can expand the modified 
Poisson equation around tpQ. If the small correction to 
the potential is denoted by tpx, one can write the expan- 



V • 



V a o 



+V • [[[f ■ Vipi]L(x )f + Vipi}fJ,(x )} + h.o. (16) 



Here L is a function given by x ^,, i f? , and xo is given by 
In the interior of the shell, the field is caused 



r 2 dr 



r 2 [l + L(x )]fi(x ) 



dr 



—C ipi 



]W1 

a 



■vih. 



for r > R 

for r < R 
(17) 



where C 2 is the angular-momentum-squared operator. At 
r = R, the jump conditions for pi should be satisfied: 



[1 + L{xq))h{xq) d ^ out lr ^ J; 



dipt 



Or 

l in 

dr 



(18) 



09 lr ~ R 89 
Let us now consider a particular form of the perturbation 



r=R = 4irG<r sh Y(9) 




\r=R 



pi = -ea sh Pi[cos(6)]8(r - R) 



(19) 



where Pi is the Legendre Polynomial of 1th order. The 
second of equation (JT7J) belongs to a class of equations 
known as p-Laplacian for p = 3 and N = 3, where N 
is the number of Euclidian dimensions. (For p-Laplace 
equation, see for example [Hj]. For the special case of 
p=N, see [Hj])- Little is known about its solution. How- 
ever, if the angular dependence of ip± is proportional to 
cos (8), by observation one can find that rcos{6) is a so- 
lution of equation (fT7|) inside the shell. Then for the par- 
ticular form of the perturbation with 1 = 1, the matching 
of the solution can be performed for inside and outside 
of the shell. 



* sh ecos(e) U_ + JL(£|) 2 r 



for r > R 
for r < R 

(20) 

where ^ sh = y/GM sh a a . 

Figure [2] shows the field |gM | computed from equation 
(|20| . the field \g s \ computed from equation (fT0|) . and the 
MOND field computed numerically. The numerical solu- 
tion and the perturbative solution match well at small e. 
The dashed line indicates the ratio of |g s | to the numer- 
icaly computed |gA/|. It shows that the ratio diverges as 
the perturbation parameter e approaches to zero. This 
is indicative of the fact that g s ~ s ^ 75 ~ whereas 
gM ~ e inside the shell. As a function of e, in MOND 
theory the perturbed field vanishes much more quickly 
than the naive expectation. 
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(a) (b) 

FIG. 2: (a) This figure plots the field at a particular point (r = 0.05i? and 9 = 45°) inside a shell with an I = 1 perturbation 
of the form given by equation (|19l as a function of the strength e of the perturbation. The strength of the numerical MOND 
field gM, the perturbative MOND field and the naive MOND field |g s | are shown. The dashed line indicates the ratio of |g s | 
to |gAf| with the scale given on the right axis. (b)This figure shows the field configuration inside the shell for a perturbation 
of order / = 1 with e = 0.1. Green arrows indicate the numerical MOND field, and red arrows indicate the difference between 
the numerically calculated MOND field and the naive MOND field, gM — gs- The sizes of the arrows are magnified by 10 for 
clarity. 




(a) (b) 

FIG. 3: (a) This figure plots the field at a particular point (r = 0.05i? and 9 = 45°) inside a shell with an I = 2 perturbation 
of the form given by equation (|19l as a function of the strength e of the perturbation. The strength of the numerical MOND 
field gA/, the perturbative MOND field and the naive MOND field |g s | are shown. The dashed line indicates the ratio of |g s 
to |gA/| with the scale given on the right axis. (b)This figure shows field configuration within the sphere for I = 2 at e = 0.1. 
Green arrows indicate the numerical MOND field, and red arrows indicate the difference of numerically calculated MOND field 
and the naive MOND field, gA/ — gs. The size of the arrows are magnified by 10 for clarity 

For the higher order perturbation of I > 2, the analytic FIG. EJb) and FIG. BJb) show the field configuration in- 
solution of even the 'simple' Possion equation (fl~7|) is not side the shell for / = 2 and I = 3 respectively. Green 
known, and we rely entirely on the numerical simulation. arrows indicate the MOND field gM, and red arrows in- 
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FIG. 4: ((a)This figure plots the field at a particular point (r = 0.05i? and 9 — 45°) inside a shell with an I = 3 perturbation 
of the form given by equation (|19l as a function of the strength e of the perturbation. The strength of the numerical MOND 
field gM, the perturbative MOND field and the naive MOND field |g s | are shown. The dashed line indicates the ratio of |g s | 
to |gA/| with the scale given on the right axis. (b)This figure shows field configuration within the sphere for I = 3 at e = 0.1. 
Green arrows indicate the numerical MOND field, and red arrows indicate the difference of numerically calculated MOND field 
and the naive MOND field, gM — g s . The size of the arrows are magnified by 10 for clarity 



dicatc the difference gM — g s . From these figures, it is 
apparent that gM << gs inside the shell. Similar to the 
/ = 1 case, the simple rescaling formula overestimates the 
field inside. FIG.^a) and FIG. Ufa) shows |g M | and |g s | 
at a single point (r=0.05R and 9 — 45°) as a function of 
perturbation parameter e. As shown in the figures, just 
like I = 1 case, the ratio |g s |/|gAf| diverges as e — * 0. It 
indicates that gM vanishes faster than g s in the limit of 
small e. 

In general, the MOND field due to the small pertur- 
bation is suppressed compared to g s , and the aspherical 
perturbation external to the shell may not influence the 
internal system as much as one might naively expect. 



SHIELDING AND ANTI-SHIELDING EFFECTS 

Since the MOND equations are non-linear, it might be 
expected that changing the mass of the shell while keep- 
ing the perturbation term constant would affect the field 
inside the shell. (In Newtonian dynamics, it does not.) 
This is indeed the case. Again, we consider a pertur- 
bation of the form equation p^|) with 1 = 1. Keeping 
the a s h = M s /R 2 term in equation (fT9)l constant while 
increasing the mass of the shell from M s to M new , the 
potential inside the shell becomes 



Iplnew — V GM n . 



ewCLQ— eCOS(fl) 



Mne 



r 6 / r t new \ 2 r_ 
2R 56 \ R ) R 



-ecos(#) 



r 

2R 



_6_ 

56 



/ ^tnew \ * 

\ R J 



R 



(21) 



Comparing equation (I2ip and equation (|20|) . equation 

([2"T) is suppressed by the factor of yj ^ Ish . Hence, the 

field within the shell is reduced by adding more mass 
to the shell. The perturbed field within the shell is re- 
duced by the spherically symmetric part of mass distri- 
bution. This can be thought as a screening effect of the 
spherical shell. As the spherically symmetric part of the 
mass distribution becomes more dominant, the perturbed 
field within the shell is reduced. The field corresponds to 
equation (f2Tj) . and is represented by the curve for R2 = R 
in FIG. [3 M new can be decomposed into M s h+M 2 where 
M2 is the spherically symmetric component that has been 
added. 

The situation described in equation (|21[) can be 
thought as a special case of a two-shell system where 
the radius of the first shell R and that of the second 
shell R2 coincide. FIG. [5] depicts the general case where 
the two radii are different. Five curves are plotted cor- 
responding to different radii for the second shell. The 
horizontal dashed line indicate the MOND field of the 
first shell with the perturbation when the second shell is 
not present. In the limit where M2 vanishes, the field 
will reduce to this value. From FIG. El one can observe 
that, for a generic value of the inner mass (M2 ~ M s h), 
the field is screened regardless of the radius of the second 
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Log(M 8 /M total ) 

FIG. 5: This figure shows the magnitude of field at r = 0.05i? 
and 9 ~ 45° as a function of M-z/MtotaU The mass and the 
radius of the first shell is fixed at R = 0.64Mpc and M s j = 
2.0 x 1O 13 M . The perturbation given in equation (fH))l with 
I — 1 and e = 0.07 is put on the first shell. Each curve 
corresponds to the second shell having a different radius. The 
horizontal dashed line indicates the MOND field of the single 
shell system with equivalent mass, radius and perturbation 
as the first shell. The scale on the right is normalized to this 
value. The dotted line represents the field corresponds to the 
perturbative solution of equation (f2"Tj) . 



shell. Increasing the mass of the spherically symmetric 
part of distribution, normally screens the field from the 
aspherical part of the mass distribution. The screening 
will be most efficient for the special case when the two 
shells coincide (i.e. R = R2), and in this case field in- 
side is always suppressed. When R2 > R, the field is 
again always suppressed. For a given value of M2, the 
suppression factor is less than that of the special case 
R = i?2- The case R2 < R exhibits a carious feature. 
The field inside the second shell can be enhanced when 
M2 << M s h, this is an anti-screening effect. We also 
note that this enhancement is larger for R2 « R. 

When M 2 « M sh and R 2 « R, the second shell be- 
longs to the local system while the first shell is the back- 
ground mass distribution. The local shell then picks up 
and enhance the perturbations of the background mass 
distribution. 

Finally, we consider the two-shell system with a per- 
turbation other than pure 1 = 1. Specifically, a spherical 
cap is removed from the first shell with opening angle a, 
measured from positive z-axis, leaving behind a spherical 
bowl. The mass density of the bowl is increased so as 
to keep the total mass of the bowl equal to that of the 
original full shell. We consider opening angles of 8.2° and 
18.6°. A perturbation of this form contains the sum of 




L°g(M ln /M lolal ) 

FIG. 6: This figure shows the magnitude of the field at r = 
0.05-R and 6 ~ 45° as a function of M%. The radius of the 
second shell is fixed at R2 = 0.3R. The perturbation on the 
outer shell is provided by an opening on top of the shell. Two 
curves are shown for the opening angle of 18.6° (blue) and 
8.2° (green). The dashed line indicates the corresponding 
field without the presence of the inner shell. 

Legendre Polynomials up to infinite oreder. 

Although the form of the perturbation is very differ- 
ent, FIG1H] shows much the same feature as FIGG3 For 
M 2 r~j M s h the field is screened. For larger breaking 
of the spherical symmetry (larger a), the screening be- 
comes less effecient and it requires more mass for the 
second shell to suppress the field. For M 2 « M s h, 
the field within the second shell is again enhanced. It 
might be expected that the enhancement becomes promi- 
nent when the generic MOND field of the second shell 
~ \/GM2~ao/R 2 becomes comparable to the field pro- 
duced by the bowl-only system. For example, for the 
opening angle of a = 8.2°, the bowl produces an aver- 
age MOND field of - 0.0052a at r = R 2 . The generic 
MOND field of the second shell becomes comparable 
when Log[M 2 / Mtotai] ~ —2.1. Similarly, for the opening 
angle of a = 18.6°, the MOND field of the second shell 
becomes comparable at Log[M2/M tota i\ ~ —1.6. These 
values reasonably estimate the mass when the enhance- 
ment starts to become important. 

CONCLUDING REMARKS 

We have shown, by examining a simple toy system of 
a shell with a perturbation that the naive MOND equa- 
tion (flO|) overestimates the field inside the shell. For a 
small perturbation on the shell signified by the pertur- 
bation paramter e, the field within the shell shell is of 
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order e and not as predicted by the naive MOND 
equation. To properly estimate the field within the shell, 
the jump condition on the surface of the shell has to be 
taken into account. We have also shown that the field 
due to the perturbation within the shell can be screened 
by the spherical distribution of the mass. The screening 
is strongest if the spherical distribution is positioned near 
the shell. This contrasts to Newtonian gravity where ad- 
dition of the spherical distribution docs not affect the 
field configuration within the distribution. Finally, the 
anti-screening effect is noticed when the small and light 
shell is added to the system. The enhancement inside 
the small shell becomes significant when the generic field 
of the shell ~ y/GMao/R becomes comparable to the 
perturbation field near the shell. The enhancement is 
generally larger for the shell with smaller radius. 

Let now suppose that, in the simplest picture, a cluster 
is modeled by a spherical shell and the external field is 
produced by an aspherical perturbation on the shell. In 
order to estimate the magnitude of the external field, 
one needs to know not only the size of perturbation but 
also the total mass of the shell. Applying naive MOND 
equation on this system to estimate the external field 
would be quite inadequate. If the galaxy can be modeled 
by a small shell within a large, more dominant shell (a 
cluster) , then the external field within the galaxy is likely 
to be enhanced. The enhancement depends on the ratio 
Mg a i/M c i u and the size of the galaxy, but as shown in 
FIG. 03 it can be of the same size as the external field 
itself. 

ACKNOWLEDGEMENTS 



group at CWRU. 



[1] De-Chang Dai, Reijiro Matsuo and Glenn Starkman, 

larXiv:0811.1565l 
[2] Mordehai Milgrom, Astrophys. J. 270: 365-370, 1983 
[3] Mordehai Milgrom, Astrophys. J. 270: 371-383, 1983 
[4] Mordehai Milgrom, Astrophys. J. 270: 384-389, 1983 
[5] Rafael Brada and Mordehai Milgrom, 

Mon.Not.Roy.Astron.Soc.276: 453-459, 1995 
[6] Kent S.M. AJ, 93:816, 1987 

[7] G.W. Angus, B. Famaey and DA. Buote 

Mon.Not.Roy.Astron.Soc. 387: 1470-1480, 2008 
[8] Benoit Famaey and James Binney, 

Mon.Not.Roy.Astron.Soc.363: 603-608, 2005 
[9] K.Begeman,A.Broeils and R.Sanders, 

Mon.Not.Roy.Astron.Soc. 249: 523-537, 1991 
[10] Jacob Bekenstein and Mordehai Milgrom, Astrophys. J. 

286: 7-14, 1984 
[11] Mordehai Milgrom, Astrophys. J. 302: 617-625, 1986 
[12] Mordehai Milgrom, Astrophys. J. 333: 689, 1988 
[13] Peter Lindeqvist, Notes on the p-Laplace equation, Re- 
port, University of Jyvskyl, Department of Mathematics 

and Statistics, Jyvskyl, 2006 
[14] Rouba Borghol and Laurent Veron, Communications in 

Partial Differential Equations, Volume 32, Issue 6:1001- 

1015, 2007 |arXiv:matn/0610016l 
[15] X.Wu, H.S.Zhao, B.Famaey, G. Gentile, O.Tiret, 

F. Combes, G.W.Angus, A.C.Robin, Astrophys. J. 665: 

L101-L104. 2007 
[16] B.Famaey, J.P.Bruneton, H.S.Zhao, 

Mon.Not.Roy.Astron.Soc.377:L79-L82, 2007 
[17] X.Wu, B.Famaey, G.Gentile, H.Perets, H.S.Zhao 

Mon.Not.Roy.Astron.Soc. 386, Issue 4, 2199-2208 



GDS and RM were supported by a grant from the US 
Department of Energy to the particle astrophysics theory 



